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ORIGINAL ARTICLE 

The spatial distribution of overweight and obesity among a birth 
cohort of young adult Filipinos (Cebu Philippines, 2005): 
an application of the KuUdorff spatial scan statistic 

DL Dahly\ P Gordon-Larsen2'^ M Emch^'^ J Borja^ and LS Adair^'^ 

OBJECTIVES: The objectives of the study were to test for spatial clustering of obesity in a cohort of young adults in the Philippines, 
to estimate the locations of any clusters, and to relate these to neighborhood-level urbanicity and individual-level socioeconomic 
status (SES). 

SUBJECTS: Data are from a birth cohort of young adult (mean age 22 years) Filipino males (n = 988) and females (n = 820) enrolled 
in the Cebu Longitudinal Health and Nutrition Survey. 

METHODS: We used the Kulldorff spatial scan statistic to detect clusters associated with unusually low or high prevalences of 
overweight or obesity (defined using body mass index, waist circumference and body fat percentage). Cluster locations were 
compared to neighborhood-level urbanicity, which was measured with a previously validated scale. Individual-level SES was 
adjusted for using a principal components analysis of household assets. 

RESULTS: High-prevalence clusters were typically centered in urban areas, but often extended into peri-urban and even rural areas. 
There were also differences in clustering by both sex and the measure of obesity used. Evidence of clustering in males, but not 
females, was much weaker after adjustment for SES. 
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INTRODUCTION 

Urbanization, modernization and globalization have radically 
altered human environments, resulting in sweeping changes to 
the way we work, play and eat.^'^ The most dramatic changes are 
taking place in the developing world, where the pace of 
urbanization has been especially rapid.^ In the wake of these 
environmental changes, obesity has emerged as a global public 
health problem,^"^ even in contexts where underweight is still 
prevalent.^ 

While it is broadly accepted that physical and social environments 
can promote obesogenic behaviors and limit healthy options,^"^^ 
more research on liow they do so is needed to inform intervention 
efforts.^ ^'^^ Responding to this need, researchers have spent the last 
two decades blending theoretical and methodological perspectives 
from various disciplines to identify environmental determinants of 
obesity (for example,^ '''^^). However, this growing area of research 
has almost exclusively focused on high-income countries, while 
related research in lower-income countries mostly consists of 
comparisons of urban vs rural obesity prevalences. 

With this gap in mind, we employed the Kulldorff spatial scan 
statistic^^ to detect clusters with unusually high or low 
prevalences of obesity among young adults living in a large 
metropolitan area in the Philippines. We then compared the 
locations of these spatial clusters with the urbanicity of the area, 
using a continuous scale measure of urbanicity that captures 
environmental heterogeneity between and within 'urban' and 



'rural' areas. We also evaluated the degree to which any clusters 
were explained by the spatial distribution of individual-level 
socioeconomic status (SES). 



MATERIALS AND METHODS 

Study site 

Metro Cebu (pop 1.9 million), on the east coast of Cebu Island in the 
central Philippines, is composed of three cities and seven municipalities. 
It is further divided into 270 administratively defined neighborhoods called 
barangays (average area 2.65 km^) comprising a 720 km^ contiguous area. 
The study area includes densely populated urban centers, less dense peri- 
urban areas, rural towns and more isolated mountain and island areas. 

Study design and sample 

Data are from the Cebu Longitudinal Health and Nutrition Survey (CLHNS), 
a community based, one-year birth cohort.^^ A single-stage cluster 
sampling procedure was used to randomly select 33 barangays. 
Pregnant women who gave birth between May 1, 1983 and April 30, 
1984 were recruited from these barangays into the study. More than 95% 
of identified women agreed to participate. A baseline interview was 
conducted among 3327 women during their sixth or seventh month of 
pregnancy. Another survey took place immediately after birth. There were 
3080 non-twin live births which make-up the CLHNS birth cohort. 
Subsequent surveys were conducted bi-monthly to age two years, then 
in 1991, 1994, 1998, 2002 and 2005. 

We used birth cohort data on barangay of residence and multiple 
adiposity measures collected in 2005 when the study participants were 
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young adults (20-22 years of age). Women pregnant in 2005 were 
excluded (n = 73). Anyone with missing data for barangay residence or any 
outcome measures were also dropped (n = 7), resulting in a sample of 988 
males and 820 females. By 2005, participants were living in 161 different 
barangays, though 77% of males and 76% of females were still located in 
the original 33 sample barangays (Supplementary Figure 1). Intra-barangay 
sample sizes were highly variable, ranging from 5 to 163 participants 
(mean: 43.8, s.d. 35.8) in the original 33 barangays; and sparse in the 
additional barangays, ranging from 1 to 15 (mean: 3.4 s.d. 2.8). 

Barangay-level data were collected for each round of the survey. The 
surveys included information on physical characteristics, infrastructure and 
utilities, social services, community organizations, industrial and commer- 
cial establishments, labor markets and wage rates. Data were obtained 
from barangay officials or other key informants. Population sizes were 
taken from the year 2000 Filipino census. 

Measures 

Neighborhoods. Barangay of residence in 2005 was used to define 
participant neighborhoods. The assumptions implied are constant 
barangay-level effects across respondents in a given barangay, and that 
the barangay is a reasonable approximation of a person's 'activity space,' 
defined as the set of locations a person encounters in the course of their 
daily activities.^^'^^ These assumptions are common in health geography. 
Barangays have their own elected officials and budgets, community 
centers and so on, and Cebuanos can easily identify their barangay. Thus, 
barangay of residence is a better measure of neighborhood, one that 
encompasses social factors and not just space, than the administratively 
defined neighborhoods used in other contexts. 

Urbanicity. While most research uses the urban-rural dichotomy to 
describe urbanicity, our research uses a continuous measure that captures 
a range of variation in urbanicity across a single dimension. The starting 
point for this measure is a previously designed urbanicity scale. Its 
description, rationale, and validation are given in Dahly and Adair.^° 

Briefly, the scale is made up of seven components derived from data 
collected for the CLHNS barangay-level surveys. The components are 
population size; population density; communications (availability of mail, 
telephone, internet, cable TV and newspaper services); transportation 
(paved road density and public transportation services); markets (presence 
of gas stations, drug stores, grocery stores and the number of small 
commercial kiosks); educational facilities; and health services. Theoretically, 
the scale represents an underlying latent construct labeled urbanicity that 
is imperfectly reflected in each of these seven components. 

Since publishing the details of the scale's creation, we've modified it by 
making the urbanicity value for a given barangay a function of its own 
score and the scores of surrounding barangays (for example, an urban 
barangay surrounded by other urban barangays will have a higher final 
score than an urban barangay with the same initial value that is bordered 
by more rural barangays). These values were created with the ESRI ArcMap 
inverse distance weighting interpolation tool, using the default settings. 
For more detail on inverse distance weighting, please see Waller and 
Gotway.^^ A map of interpolated urbanicity values is given in Figure 1. 

Anthropometrics. We used three anthropometric measures to define 
overweight (OW) or obesity. They were body mass index (BMI), waist 
circumference (WC), and percent body fat (BF%). All anthropometrics were 
collected by trained field staff during in-home interviews using techniques 
described in Lohman et ai?^ Weight was measured with calibrated 
portable scales to the nearest g. Height was measured with a folding 
stadiometer to the nearest tenth of a cm. BMI was calculated as measured 
weight (kg) divided by measured height (m) squared. BMI is the standard 
measure of adiposity for public health research^^ because it is strongly 
correlated with both BF% and total fat mass,^^ but it does not differentiate 
fat mass from lean mass. BMI was dichotomized as OW or obese based on 
respective cut points of ^ 23 kg m ~ ^ and ^ 25 kg m ~ ^, primarily based on 
evidence that lower cut points in Asian populations may be more 
appropriate for describing health risks (for example, Misra^^ and Gallagher 
et al}'^), and better reflect an important 'public health action point'.^° 

WC was measured to the nearest mm, at the midpoint between the 
bottom of the ribs and the top of the iliac crest, using a measuring tape. 
WC is a measure of central obesity, which is thought to be an important 
driver for a number of important metabolic disorders^^ compared with 
other patterns of fat distribution. High WC (HWC) was defined as 
WC>85cm in males or >80cm in females.^^ 



BF% was calculated as described by Durnin and Wormersley^^ using the 
sum of three skin folds (triceps, subscapular and suprailiac), measured to 
the nearest mm using Lange calipers. It is a more valid measure of 
adiposity than BMI, though it is typically measured less reliably. High BF% 
(HBF) was defined as BF% ^25 in the males, or ^38 in the females.^"^ 

Socioeconomic status. To capture SES, we used principal components 
analysis to evaluate a set of household assets measures (for example television, 
land and so on), a method commonly used with Demographic and Health 
Survey data. The measure used here is the first principal component 
(eigenvalue 4.2) from this analysis, and has been used in previous analyses.^^'^^ 

Analytical methods 

A spatial cluster can be defined as a contiguous geographic space for 
which the value of some characteristic is unusual when compared with 
the space surrounding it. The characteristic we are interested in is the 
prevalence of OW or obesity in our sample. To detect clusters, we used the 
spatial scan statistic as implemented by the software SaTScan and 
employed the Bernoulli model^^'^^ for binary outcomes. 

We derived Cartesian coordinates (x,y) for the center-point of each 
barangay using ArcGIS. Individuals in the study were assigned the location 
of their respective barangay's center-point, resulting in a set of point 
locations characterized by the number of 'cases' of each outcome and the 
total number of study participants at that location. 

For each of these locations, SaTScan draws all possible contiguous scan 
windows, centered on that point, that contain ^50% of the total study 
population. For each of these windows, a prevalence (pj is calculated as the 
number of observed cases divided by the total population residing within 
the window. This is compared with the observed prevalence [q) for 
participants residing outside the given window, resulting in a prevalence 
ratio (PR = p/(gj. The null hypothesis tested for each window is Hq: p/q^ 1. 
In this analysis, we set out to detect clusters of higher and lower than 
expected prevalence, thus the alternate hypothesis is H^ p/q^^l. 

The goal is to identify the window with a prevalence ratio that is the 
least likely to occur given Hq. To do this, a likelihood is calculated for each 
window. For the Bernoulli model detecting both high and low-prevalence 
clusters, the likelihood function is: 

W \ n ) \N-n) \ N-n ) 

where c is the number of cases in the window, C is the total number of 
cases in the sample, n is the number of observations in the window, and N 
is the total number observations in the sample.^^ The likelihood function is 
maximized over all the windows, and the one with the maximum 
likelihood is identified as the cluster that is least likely to have occurred by 
chance.^^ This cluster is the 'primary cluster.' 

To obtain a P-value, the same analysis is repeated on 999 random 
replications of the data generated under the null hypothesis, and the 
maximum likelihood from the real data is ranked (/?) along with maximum 
likelihoods from each of these Monte Carlo simulations.^^ The P-value of the 
least likely cluster is given by R/1000. For this analysis we have elected to 
report all primary clusters for which Hq can be rejected at P<0.15, and any 
secondary clusters with P<0.15 that don't geographically overlap with the 
primary cluster. Highly unusual clusters, given Hq, are indicated by P<0.05. 

The shape of the scan window used to detect clusters can take on a 
variety of forms. We chose to use an elliptical scan window^^ because the 
urban core of Metro Cebu is elongated, running roughly southwest to 
northeast (see Figure 1). By hypothesizing that adiposity clusters will 
geographically coincide with the most urban areas of Metro Cebu, we are, in 
effect, hypothesizing that the clusters will be elliptical in nature. The elliptical 
scan window has slightly higher power for long narrow clusters, and slightly 
lower power for circular or more compact clusters.^^ We used a medium 
strength non-compactness penalty which favors less eccentric clusters.^^ In 
practice, this penalty helps prevent the detection of very long, thin clusters 
that are artifacts driven by high prevalences at the two ends of the ellipse. 

We then repeated our analysis, adjusting for individual-level SES. The 
goal was to evaluate the degree to which spatial clusters of adiposity are 
explained by the spatial distribution of SES among study participants. In a 
previous analysis using generalized estimating equations, we found strong 
positive associations between adiposity and multiple indicators of SES in 
the males, but not in the females.^^ The adjustment for SES was made by 
stratifying the samples by sex-specific tertiles of SES and using the multiple 
data sets option in SaTScan. More detail on how this adjustment is 
calculated can be found in the SaTScan user's guide.^^ 
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Figure 1. Urbanicity in Metro Cebu, 2005. To further validate that the urbanicity scale provides a valid description of the rural-urban gradient 
in Metro Cebu, we connpared its spatial distribution to a Landsat 7 ETM+ innage of the study area (http://landsat.gsfc.nasa.gov/), which is a 
false color connposite that depicts vegetation as shades of red, while nnixes of bare soil and innpernneable land cover (buildings, roads and so 
on) appear green. We also connpared our innage with a SRTM elevation nnap (http://srtnn.csi.cgiar.org/) because nnountains constrain urban 
developnnent to the northwest. The nnap depicts higher elevations as darker shades of blue, because nnountains constrain urban developnnent 
to the northwest. 



Clusters are displayed by outlining the barangays that make-up the 
cluster (in red for high-prevalence clusters and in blue for low-prevalence 
clusters), then combining that information with urbanicity maps of 
the study area using ArcGIS. Urbanicity maps were created using the 
continuous urbanicity scale, as well as the barangays' urban-rural 
designations from the 2000 Philippines census. 



RESULTS 

Table 1 sunnmarizes the distribution of adiposity measures in our 
sannple. The sannple is lean overall, but characterized by wide variation 
in nneasures, indicative of a population suffering fronn a double 
burden of over- and under-nutrition (consistent with Pedro et al."^^). 

Annong nnales (Table 2), there were high- and low-prevalence 
clusters detected for each outconne except HWC, for which there 
was only a high-prevalence cluster. Generally, the prevalence of a 
given outconne annong nnales residing within a high-prevalence 
prinnary cluster was nnore than twice that of nnales living outside of 
the cluster. The high-prevalence clusters for OW and obese were 
highly unusual (P<0.05) given the null hypothesis of connplete 
spatial randonnness. Conversely, the high-prevalence clusters for 
HBF and HWC were not highly unusual (P = 0.104 and 0.108, 
respectively). In every instance, these high-prevalence clusters 
were located in the urban core of Metro Cebu, with nninor 
variations in the set of barangays contained within thenn 



Table 1. Adiposity measures in 988 male and 820 female young adults 
(mean age 22 years) enrolled in the Cebu Longitudinal Health and 
Nutrition Survey, 2005 




Moles (n = 988) 


Females (n = 820) 




BMI WC BF% 
kgm~^ cm 


BMI WC BF% 
kgm~^ cm 


Mean 
s.d. 

Median 
Maximum 
Minimum 
% OW^ 
% Obese^ 
% HBF^ 
% HWC"^ 


21.03 72.17 16.70 
3.05 7.53 5.10 
20.46 70.8 16.07 
40.33 112 33.59 
14.48 56.5 3.69 
19.7 (17.3-22.3)^ 

9.4 (7.7-11.4) 

6.7 (5.2-8.4) 

6.1 (4.7-7.7) 


20.23 67.91 32.65 
3.17 7.45 4.77 
19.67 66.5 33.05 
41.17 112.2 45.56 
13.93 53.7 17.09 
15.4 (13.1-18.0) 

7.8 (6.2-9.8) 
12.2 (10.0-14.6) 
6.5 (5.0-8.4) 


Abbreviations: BF%, percent body fat ; BMI, body mass index; HBF, high % 
body fat; HWC, high waist circumference; OW, overweight; WC, waist 
circumference. ^OW defined as BMI ^ 23 kg m " ^, "^Obese defined as 
BMI^25kgm ~^. ^HBF defined as BF% ^25 in males or ^38 in females. 
'^HWC defined as WC>85cm in males or >80cm in females. 
^95% confidence intervals for proportions reported as estimate (lower 
limit to upper limit), and calculated using the Wilson procedure.^^ " 
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Table 2. SaTScan^ clusters for multiple anthropometric measures in young adult males (A/ = 
Cebu, 2005 


988) and females (A/ = 820) enrolled in the CLHNS; Metro 


Outcome 


Total cases 


Total prevalence 








Cluster information 










ID 


n 


cases 


pb 




P-value 


Males 


















Obese^ 


93 


0.09 


1 ° 


95 


0 


0.00 




0.015 








2° 


41 1 


58 


0.14 


2 33 (^ 56 to 3 47) 


0.041 


OW 


195 


0.20 


1° 


492 


134 


0.27 


1.T\ (1.68 to 2.92) 


0.001 


HWC^ 






2° 


95 


3 


0.03 


0.15 (0.05 to 0.45) 


0.002 


60 


0.06 


r 


428 


41 


0.10 


2.82 (1.66 to 4.79) 


0.104 


HBF^ 


66 


0.07 


r 


202 


2 


0.01 


0.12 (0.03 to 0.49) 


0.021 








2° 


478 


47 


0.10 


2.63 (1.57 to 4.43) 


0.108 


Females 


















Obese 


64 


0.08 


r 


405 


47 


0.12 


2.83 (1.65 to 4.84) 


0.062 


OW 


126 


0.15 


r 


407 


85 


0.21 


2.10 (1.48 to 2.97) 


0.020 


HWC 


53 


0.06 


r 


405 


42 


0.10 


3.91 (2.04 to 7.49) 


0.010 


HBF 


100 


0.12 


1° 


214 


43 


0.20 


2.13 (1.48 to 3.07) 


0.125 


Abbreviations: HBF, high % body fat; HWC, high waist circumference; ID, cluster identifier (1° primary cluster; 2° secondary cluster); OW, overweight; 
P, prevalence; PR, prevalence ratio. ^SaTScan Bernoulli model; elliptical scan window; medium non-compactness penalty; maximum cluster ^50% of the 
population; overlapping clusters not reported; P-values two sided. "^Prevalence defined as (cases/n) within a cluster. ^Ratio of the prevalence inside a cluster to 
the prevalence outside of the cluster. 95% confidence intervals reported as estimate (lower limit to upper limit). ^Obese defined as BIV1I^25 kg m"^. 
^OW defined as BIV1I^23 kg m"^. ^HWC as WC>85cm in males or >80cm in females. ^HBF defined as BF% ^25 in males or ^38 in females. 



(see Figure 2a). Low-prevalence clusters, where outcomes were 
virtually non-existent, occurred in the rural south or southwestern 
regions of Metro Cebu, though the low-prevalence cluster of HBF 
extended into some urban areas. 

Like the males, a high-prevalence cluster was detected for all 
four outcome measures in females (Table 2), though the obese 
and HBF clusters were not highly unusual given the null 
hypothesis of complete spatial randomness (P = 0.062 and 0.125, 
respectively). The prevalence of a given outcome among females 
living in the respective cluster was generally more than twice that 
of females living outside the cluster, but almost four times larger 
for the HWC cluster. Each of these high-prevalence clusters was 
located in the urban core of Metro Cebu (see Figure 2b). There 
were no low-prevalence clusters detected among the females. 

To illustrate how informative the urban-rural dichotomy would 
be for describing cluster locations, we redisplayed the male 
clusters on a map of Metro Cebu barangays defined by the urban- 
rural dichotomy (Figure 3). For example, though the magnitude of 
the prevalence ratio comparing obesity in urban versus rural 
males (PR 1.71; 95% CI 1.05 to 2.79) is similar to that found in the 
high-prevalence obesity cluster (PR 2.33; 95% CI 1 .56-3.47; top left. 
Figure 2a), if we only intervened in the urban barangays, we would 
be ignoring 18 rural barangays included in the obesity spatial 
cluster (a total population size of 96,472 according to the year 
2000 Filipino census), and acting in 30 urban barangays not 
included in the spatial cluster (total population size 380,322). 
Furthermore, looking at the HBF clusters (bottom right. Figure 2a), 
intervening to reduce obesity in urban areas would target 
barangays that are in fact part of a low-prevalence cluster. 

After adjustment for individual-level SES (Supplementary 
Table 1), there was no evidence of spatial clustering of obesity, 
HWC or HBF in the males. There was still strong evidence of a 
high-prevalence cluster of OW in the males, but no evidence for a 
low-prevalence cluster. Adjustment for SES did not appreciably 
impact the spatial clusters in the females. 



DISCUSSION 

Motivated by 



the increasing prevalence of obesity in the 



Philippines,"^''"^^ and other lower-income countries,^ particularly 
in urban areas, we aimed to identify spatial clusters of obesity in 



Metropolitan Cebu. By identifying spatial clusters and relating 
them to the urbanicity of the study area, we hoped to facilitate 
further etiological research aimed at identifying environmental 
determinants of obesity in this study area. We also hoped to more 
directly inform public health practice by explicitly identifying areas 
where obesity is most common, and by evaluating how useful 
simple urban-rural classifications are for identifying these areas 
when spatial data are not available. 

Cluster locations were consistent with the idea that urban areas 
are obesogenic. High-prevalence clusters were centered on the 
highly urban core of Metropolitan Cebu, and low-prevalence 
clusters were found in rural areas. However, high-prevalence 
clusters typically covered substantial proportions of the study area 
that included both urban and peri-urban areas. The observed scale 
of the clusters likely reflects the simple fact that barangay of 
residence is not a complete description of a person's activity 
space. The lack of smaller clusters also suggests that there are not 
any environmental determinants of obesity that exert strong, 
highly localized influence. 

Evidence of spatial clustering varied by sex and outcome 
measure. The most striking difference between the sexes was that 
evidence of clustering for HWC was very strong for females 
(P = 0.010), but not for males (P = 0.104). Furthermore, the female 
HWC cluster was characterized by a PR roughly twice that of the 
other high-prevalence clusters (PR ~4 vs ~2). These results are 
consistent with a previous analysis of the same sample in which 
we observed gender differences in intra-class correlations from 
multi-level models describing the proportion of variance in 
obesity and HWC described at the barangay level (males: 0.18 
and 0.10, respectively; females: 0.05 and 0.18).^^ These results raise 
the interesting possibility that central fat distribution and total 
adiposity each have unique environmental determinants that vary 
by gender. 

High-prevalence clusters were not detected for all of the 
outcomes considered. In the males, there were clusters for both 
obesity and OW, but not HWC or HBF. We cannot rule out that 
these clusters are reflecting influences on both lean and fat mass. 
For females, a lack of clustering in HBF and obesity lead to similar 
reservations. While BF% is a more valid measure of adiposity than 
BMI, it is measured less reliably due to its reliance on skin folds, 
versus height and weight, which are measured much more 



Nutrition & Diabetes (2013) 1 -8 



© 2013 Macmillan Publishers Limited 



Spatial distribution of overweight and obesity 
DL Dahly et al 



a N 



Urtanicity 
^■oto15 (Rural) 
I I 15+ to 20 
I I 20+ to 25 
I I 25+ to 30 
I I 30+ to 35 
I I 35+ to 40 
40+ to 45 
I I 45+ to 50 
I 1 50+ (Urban) 

Sample 
barangay 

I 1 

High 

prevalence 
cluster 

I 1 

Low 

prevalence 
cluster 




Urtanicity 
^1 0to15 (Rural) 
I I 15+ to 20 
I I 20+ to 25 
I I 25* to 30 
I I 30+ to 35 
r I 35+ to 40 
40+ to 45 
I I 45+ to 50 
I 1 50+ (Urban) 

Sample 
barangay 



High 

prevalence 
cluster 

I — I 




Figure 2. (a) High and low-prevalence obesity clusters in CLHNS males, n = 988 (Metro Cebu, Philippines, 2005). (b) High-prevalence obesity 
clusters in CLHNS females, n = 820 (Metro Cebu, Philippines, 2005). BF%, percent body fat; BMI, body mass index; WC, waist circumference. 
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Figure 3. Male obesity clusters displayed with barangays designated as urban or rural by the Philippines census (Metro Cebu, Philippines, 
2005). BF%, percent body fat; BMI, body nnass index; WC, waist circunnference. 



accurately. Thus measurement error in BF% could have reduced 
our power to detect any such clusters. WC is also measured less 
reliably, though this did not prevent the detection of a high- 
prevalence cluster in the females. 

We also cannot rule out that the existence of clusters are driven 
by compositional factors (that they disproportionately contain 
people more susceptible to obesity) rather than contextual 
causes. However, because the sample is a one-year birth cohort, 
confounding of the space-outcome relationship by age is not 
possible. Furthermore, the ethnic and cultural make-up of Cebu is 
fairly homogenous. Based on previous research on the socio- 
economic determinants of obesity in this sample, we did adjust for 
a measure of household assets. We found that evidence for 
clustering in the males was much weaker after this adjustment, 
though there was no impact of the adjustment on the female 
results. This finding was consistent with a large body of evidence 
suggesting that SES is an important determinant of obesity in 
males but not females in lower-income contexts.^^""^^ 

In addition to gaining some etiological insight, understanding 
the spatial distribution of disease can also aid public health efforts 
by providing the best possible information on where to focus 
interventions. Unfortunately, the spatial data required to detect 
spatial clusters of obesity are often not available. In lower-income 
countries, this information is typically approximated by looking at 
obesity outcomes between areas defined as 'urban' or 'rural.' As 
obesity intervention efforts in developing countries increase, it 
seems likely they will rely on these urban-rural differences to 
target interventions. 

Our results suggest that targeting obesity interventions in this 
manner should be done with caution. Though our study area is 
the second largest metropolitan area in the Philippines, there is a 



great deal of environmental heterogeneity in Metro Cebu and 
targeting interventions at the entire area could be inappropriate. 
Furthermore, targeting the administratively defined urban areas 
within Metro Cebu also seems problematic because this descriptor 
does not accurately describe where the highest prevalences of 
obesity are. This information is valuable because lower-income 
countries have fewer public health resources, but must contend 
with a greater variety of public health problems. 



Strengths and limitations 

For this analysis, we took an urban health perspective"^^ to clearly 
link our work to previous research in lower and middle income 
contexts that has focused on crude urban-rural comparisons of 
obesity prevalences (for example, Mendez et al7). These studies 
have typically found that urbanites are more likely to be obese 
than their rural counterparts, though the degree of difference 
varies widely between studies (for example, Mendez et al7). 
However, these studies often compare urban and rural areas that 
are located in completely different regions of a country; there are 
no consistently used definitions of urban and rural, and the use of 
the urban-rural dichotomy to describe environments likely 
obscures a great deal of environmental heterogeneity.^^ To help 
overcome some of these limitations, we employed a previously 
developed continuous measure of urbanicity that better described 
the environmental heterogeneity in Metro Cebu. 

While studies of the environmental determinants of obesity are 
becoming more common, most rely on multi-level models (for 
example, Rundle et al^^) that do not account for spatial 
dependencies at the higher level, or explicitly link outcomes to 
space (but rather to characteristics of space). This approach 
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requires strong a priori hypotheses about which specific 
environmental characteristics you think are important, while our 
analysis was more exploratory. Given our goals, an explicitly 
spatial cluster analysis was more appropriate. 

While there were many available cluster detection methods, the 
Kulldorff spatial scan statistic seemed ideal in that it both locates 
clusters and provides a statistical test of how unusual the clusters 
are given the null hypothesis of complete spatial randomness. 
Furthermore, the method accounts for multiple testing; it does not 
require a priori decisions regarding the scale of the clusters; the 
clusters are robust to the spatial distribution of cases and controls 
within them; it allows for covariate adjustment; and it is 
implemented with freely available software. Lastly, because 
clusters are based on the distribution of individual study 
participants, rather than barangay-level summary measures (for 
example, proportions), participants in sparsely populated baran- 
gays do not disproportionately contribute to the clusters. 
Additional methodological comparisons and references can be 
found in Kulldorff.^^ Despite these strengths, we found few 
examples where the spatial scan statistic was used to identify 
obesity clusters (for example, Huang et al!^^ and Edwards et al. ), 
and none that were relevant to the scale and context of our 
research. 

This was the first survey in this cohort that we were able to 
assess adult obesity. Consequently, this analysis was limited by 
its cross-sectional nature. We were unable to evaluate the degree 
to which time spent in an urban area as an adult impacts obesity 
risk (such as Sobngwi et al.^^). We tried to account for the 
possibility of residential selection by limiting our sample to 
individuals who still lived in the barangay they were born in. We 
found that this exclusion had no impact on our results (not 
reported). Another potential limitation of the study is that there 
are relatively few obese people in this sample of young adults. 
However, the mothers of this birth cohort underwent a 
remarkable sixfold increase in OW and obesity prevalence 
between 1983 and 2005."^^ Space-time clustering of obesity in 
the cohort mothers will provide additional insights and help 
overcome these limitations. 



CONCLUSIONS 

Our analysis found evidence for spatial clusters of OW and obesity 
that varied by sex and outcome measure. Generally, high- 
prevalence clusters were centered on the urban core of 
Metropolitan Cebu, but extended into peri-urban areas. Low- 
prevalence clusters, where outcomes were virtually non-existent, 
were largely restricted to rural areas. The most striking result was 
the existence of a high-prevalence cluster of HWC in the females, 
characterized by a prevalence ratio of ~4. This was contrasted by 
no evidence of clustering of HWC in the males. Sex differences 
were further highlighted by the result that adjustment for 
household assets explained spatial clustering of most outcomes 
among the males but not the females. Lastly, our results suggest 
that using simple urban-rural classifications to identify areas of 
greatest risk can possibly lead to substantial amounts of 
misclassification. 
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